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Рассматривается модель стохастической волатильности. Параметры модели оцениваются методом последовательного анализа. 
Осуществляется моделирование котировок финансового инструмента на примере европейской валюты. 


Введение 

В настоящее время одним из самых актуальных 
направлений прикладной математики является по- 
строение математических моделей, адекватно описы- 
вающих эволюцию таких финансовых инструментов 
как акции, облигации, опционы, котировки валют и 
т. д. В связи с этим становится необходимым изуче- 
ние свойств, вычисление параметров и определение 
вида распределения некоторого стохастического про- 
цесса, лежащего в основе рыночных флуктуаций. 

Известно, что у эмпирической функции плот- 
ности распределения, построенной на основе исто- 
рических данных, существует ненулевой эксцесс и 
асимметрия. Кроме того, присутствует вытянутость 
функции плотности в е-окрестности точки матема- 
тического ожидания, а так же наблюдаются так на- 
зываемые «толстые хвосты», когда вероятность 
значительных изменений ценовых приращений 
выше, чем для нормального распределения. 

Для изучения распределений таких случайных 
процессов и прогнозирования будущего уровня их 
волатильности (изменчивости) было предложено 
большое количество моделей. Наиболее изученной 
является модель Блэка-Шоулса для опционов. Дан- 
ная модель позволяет построить распределение при- 
ращений временного ряда, которое уже будет иметь 
более толстые хвосты по сравнению с нормальным 
распределением и, как следствие, дает возможность 
рассчитать волатильность более точно. Тем не менее, 
хвосты остаются достаточно «тонкими» по сравне- 
нию с распределениями данных реальных финансо- 
вых рынков, что недопустимо. 

В связи с вышесказанным в настоящее время 
получили развитие другие подходы. Наиболее ин- 
тересными из них являются методы с использова- 
нием фрактального броуновского движения, когда 
временной ряд представляется в виде фрактала (и, 
как следствие, становится известен закон постро- 
ения с регулируемыми величинами начальных и 
центральных моментов); алгоритмы подгонки 
функции распределения под распределение значе- 
ний временного ряда, а так же процесс прогнози- 
рования волатильности с помощью негауссовских 
моделей типа АКСН (модели с авторегрессионной 
условной неоднородностью) [1-4), позволяющих 
проводить анализ кореллированных и высокоча- 
стотных данных. 


В данной работе рассматривается модель сто- 
хастической волатильности (, 8іосНа$ііс Ѵоіаіііііу - 
8Ѵ) первого порядка 8Ѵ(1). Для оценивания пара- 
метров этой модели в литературе предлагаются ме- 
тод максимального правдоподобия и обобщенный 
метод моментов [см., например, 1]. В работе [5] 
приводится достаточно подробное описание и 
сравнение этих двух методов применительно к 
8Ѵ-моделям. В данной работе предлагается постро- 
ить оценки неизвестных параметров методом по- 
следовательного анализа. Преимущество данного 
метода состоит в том, что оценки параметров мож- 
но получить с априорно заданной точностью. 


Постановка задачи 

Пусть 


- значение цены акций или обменного курса валют 
в момент времени я, тогда 

Н„= 1п§4 

В классе рассматриваемых моделей величина 
нормированных цен Н„ определяется как 

Н„ = у 1 +... + у„, 
где 



Модель стохастической волатильности опреде- 
ляется как [1]: 

Уп = а п 8 п ’ 

где у„ - последовательность наблюдений, е„ — стан- 
дартная гауссовская последовательность, <т„ — сто- 
хастическая волатильность. Предполагается, что 

ст; = ехр( Н „ ), Н„ = а 0 + аД_,+ ѵ„ . 

Здесь ѵ„ - стандартная гауссовская последова- 
тельность, которая предполагается независимой от е„. 

Тогда, процесс И„ будет стационарным при усло- 
вии |а,|<1, причем 
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Очевидно, что 
У 2л 


х 4 е ' 2 Эх = 


У 2л 


\х ъ е* гі1 й(х 1 /2) = 
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Еу 2 п =Еа 2 Ее 2 =Ео 2 , 

Еу л п = Е(у 4 Ее 4 = ЗЯо-;. 


Тогда коэффициент эксцесса имеет вид: 
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где Су 1 - это обратная матрица к (^(которая по 
предположению существует). Оценка Юла-Уокера 
используется здесь вместо оценки, полученной ме- 
тодом наименьших квадратов, так как она является 
асимптотически несмещенной в предположении, 
что шум ?7„^0 в наблюдениях (1). Процедура после- 
довательного оценивания в данном случае прово- 
дится в два этапа. 

Этап 1. 

Выберем такую неубывающую последователь- 
ность с к , что с*>1 и ІіітіСд=+ос, при этом 

х- 1 

Х- <о °. 

к>\ Щ 


Положительность величины К свидетельствует 
о том, что плотность распределения величин у„ в 
окрестности среднего значения вытянута вверх 
сильнее и имеет более тяжелые хвосты, чем соот- 
ветствующее нормальное распределение. 

Очевидно, что 

1п у 2 = к„ + 1п е 2 , 

Известно, что [2]: 

Е( 1о § г;) = -1,27; АЮ§г;)=4,93. 


По мере наблюдения процесса {у„}, последова- 
тельно находятся марковские моменты т к =т(с к ) по 
формуле [4]: 

т(к) = іпі{т > 2 : т_ 2 |р > к > 0}. 

2 

Здесь || I - евклидова норма для векторов и ма- 
триц, т. е. ||5|р=1г5-5 т . Для каждого т к , к> 1 модифи- 
цированная оценка Юла-Уокера определяется как: 

Ѳ(т(к)) = д;І к) [ ^ г п _ 2 - 2 п + ф)-2 ику2 - 2 ^ к) ], (4) 

п-2 


Общие положения 

Рассмотрим возможность построения оценки 
методом последовательного анализа [6, 7]. Модель 
8Ѵ может быть приведена к системе линейных ура- 
внений без потери информации: 

г „=*„+Д,> (1) 

х„=а 0 + а,х„_і+ѵ„, ( 2 ) 


где 

х п = 1о ё,К, 2 „ = 1°8У, 2 - Р, Ч, = !о§г„ 2 - Д. 

Если шумы ѵ,„ в„ являются гауссовскими, т. е. 
ѵ„,е„~У(0,1), тоД=-1,27. 

Возвращаясь к (1), (2), имеем: 

г п =а о + а і г »-і + <йі ’ ( 3 ) 

где 

<И„ =Ч„-а І Ч„- І +ѵ я . 

Это уравнение может быть переписано в виде: 

*,=2 я Г -,-Ѳ + $., 

где 

=0,г„) т , Ѳ=[а 0 ,а І ] т . 

Гарантированная оценка Ѳ для строится на ос- 
нове оценки Юла-Уокера, заданной следующими 
формулами: 

Ѳ =С-'У 2 ,2 , 

N N / і п-2 п ’ 
п-2 


г(Л)-1 

Ум = [ X ^п-2 ' ^п-\ + Ч-О 1 ) ' (И у! ' » )-і1> 

п-2 

где /и(к) - корректирующий множитель, един- 
ственным образом определяющийся уравнением: 

ДЙ 4 ,, ||2 II м2 

||^н-2 1 + Р(^)'|-^ Г (/,)_2| -Е 

п-2 


Этап 2. 


Для того чтобы закончить процедуру, введем па- 
раметр Н> 0, характеризующий точность оценива- 
ния, и определяющий последовательный план 
(ЩН),Ѳ\Н)), который состоит из длительности 
ЩН) и оценки Ѳ’(Н): 

ЩН) = 7(0+1, 


Ѳ\Н) 


а 

& 

К 7=1 У 
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Здесь ст - количество оценок типа (4), вошед- 
ших во взвешенное среднее: 


а =а(Н ) = іпЕ 



где ^=Ь(с) определяется следующим образом: 


Ь(к) 


/Г 2 



О, 


^Уь) >0, 

АсХ(} г(Н ) = 0. 


В работе [6] было показано, что для любого Н > О 
последовательный план ЩН),Ѳ’(Н ) обладает сле- 
дующими свойствами: 
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Рисунок. Сравнение исторических и модельных значений котировок курса евро (дневные цены закрытия) за период с 1 июня 
2004 г. по 21 апреля 2006 г. 


1. ЩН )< со с вероятностью почти наверное; 

2. зирЕ ||Ѳ*(//) - ѳ|| 2 < р = 2(а: + 2а 2 п )Х-. 

ТТ к> 1 с к 


Анализ эмпирических данных 

Применим метод стохастической волатильно- 
сти для моделирования временного ряда у„ котиро- 
вок европейской валюты, рисунок. 

Отметим, чтоу,=Іп(5„)— Іи(6' я _|), л= 1 ,2,3..., обла- 
дают следующими параметрами: среднее - 
1,41Ы0' 4 , дисперсия - 4,32- ІО' 5 , коэффициент 
асимметрии - 0,092, куртозис (эксцесс) - 1,747. 

Параметры модели для данного временного ряда: 
ао=0, 000561 и а, =0,998712; похожие результаты для 
котировок японской йены были получены в работе [ 1 ] . 

После моделирования прогнозных данных с ис- 
пользованием модели стохастической волатильности 
была рассчитана относительная погрешность 5„ меж- 
ду исходными данными и полученным прогнозом: 
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<5 п = 5^ 5 л - 1 . 

величина 5„ не превосходила 2,5 %. 

Анализ рисунка дает основание сделать вывод о 
том, что модель стохастической волатильности 
удовлетворительно описывает исходные данные и 
позволяет моделировать их значения с небольшой 
погрешностью. 

Выводы 

Получены оценки параметров в модели стохас- 
тической волатильности методом последователь- 
ного анализа на основе оценки Юла-Уокера. Важ- 
ными свойствами этой оценки является то, что она 
строится с априорно заданной среднеквадратиче- 
ской точностью и процедура построения оценки 
сходится при любых начальных значениях параме- 
тров модели. Кроме того, 8 V- модель использована 
для имитационного моделирования котировок 
курса евро. Показана высокая эффективность мо- 
дели стохастической волатильности. 
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